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Reflection and noise in the low spectral state of GX 339-4 
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Abstract. We analyze RXTE/PCA observations of GX339-4 in the low spectral state from 1996-1997 and show 
that the pattern of its spectral and temporal variability is nearly identical to that of Cyg X-1. In particular, a 
tight correlation exists between the QPO centroid frequency and the spectral parameters. An increase of the QPO 
centroid frequency is accompanied with an increase of the amplitude of the reflected component and a steepening 
the slope of the underlying power law. Fourier frequency resolved spectral analysis showed, that the variability 
of the reflected component at frequencies higher than ~ 1 — 10 Hz is suppressed in comparison with that of the 
primary emission. 
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1. Introduction 

It was found recently for a large sample of Seyfert AGNs 
and several observations of the Galactic X-ray binaries 
that the amplitude of the reflected component is gener- 
ally correlated with the slope of the primary power law 
emission ( ^dziarski ct al. 1999 ). Based on the numerous 
RXTE/PCA observations of Cyg X-1 [Gilfanov et al. 1999 
(hereafter Paper I) showed that this correlation is strong 
for multiple observations of this source and that the spec- 
tral parameters are also tightly correlated with the char- 
acteristic noise frequency. In particular an increase of the 
QPO centroid frequency is accompanied with a steepening 
of the slope of the Comptonized radiation and an increase 
of the amplitude of the reflected component. Studying fast 



variability of the reflected emission Revnivtsev et al. 1999 



showed that its amplitude is suppressed with respect to 
that of the primary emission at the frequencies higher than 
- 1 - 10 Hz. 

a bright and well studied X-ray 
usually classified as a black hole 
in many aspects is very simi- 
(see e.g. Tanaka&Lewin 1995 



GX339-4 is 
binary. It is 
candidate and 
lar to Cyg 



X-1 



Trudolyubov et al. 1998"i 



Zdziarski et al. 1998 



Wilms et al. 1999i [Nowak et al. 1999| ). The investi- 
gations of the connections between the spectral and 
timing properties o f Cyg X-1 (e.g. Gilfanov et al. 1999 ) 
and GX 339-4 (e.g. [Ueda et al 1994^ indicate that these 
sources could be similar from this point of view also. 
In this paper we expand the analysis of correlations 
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between spectral and temporal characteristics of the 
X-ray emission of GX 339-4 in the low spectral state 
using the Rossi X-ray Timing Explorer data and show 
that this source demonstrates the same behavior that was 
previously observed from Cyg X-1. 

2. Observations, data reduction and analysis 

We used the publicly available data of GX 339-4 observa- 
tions with RXTE/PCA from 1996-1997 performed during 
the low spectral state of the source. Our sample includes 
23 observations from the proposals 10068, 20056, 20181 
and 20183 with a total exposure time of ^^130 ksec (Table 
^). Only observations from the proposal P20183 had suffi- 
cient energy and timing resolution to perform Fourier fre- 
quency resolved spectral analysis. Therefore the frequency 
resolved analysis was carried out only for ^ 61 ksec of the 
data. 

The data screening was performed following the RXTE 
GOF recommendations: offset angle < 0.02°, Earth el- 
evation angle > 10°, electron contamination value (the 
"electron ratio") for any of PCUs < 0.1. The data from 
all PCUs were used for the analysis. The energy spec- 
tra were extracted from the PCA mode "Standard 2" 
(128 channels, 16 sec time resolution) and averaged over 
each observation. Fourier frequency resolved spectral anal- 
ysis used "Good Xenon" data (256 energy channels, 1/zs 
time resolution). The response matrixes were built using 
standard RXTE FTOOLS 4.2 tasks (Jahoda 1999). The 
background spectra for the conventional spectral analy- 
sis were constructed with the help of the "VLE" based 
model (Stark 1999). The background contribution to the 
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Fig. 1. Typical power density spectrum of GX 339-4 in 
the low spectral state. The QPO peak at the frequency of 
~0.35 Hz is shown by the arrow. In order to improve the 
statistics the power density spectrum was averaged over 
all observations of the proposal P20183 having somewhat 
different QPO centroid frequency, resulting in a somewhat 
broader profile of the QPO peak. 



frequency resolved spectra is negligible in the frequency 
and energy ranges of interest. A uniform systematic un- 
certainty of 0.5% was added quadratically to the statis- 
tical error in each energy channel. The value of system- 
atic uncertainty was chosen basing on the deviations of 
the PGA Crab spectra from a power law model (see e.g. 
Wilms et al. 19991 ). 

The energy spectra were fit in the 3-20 keV en- 
ergy range with a spectral model identical to that of 
Paper I. The model consisted of a power law without 
high energy cutoff with superposed continuum, reflected 
from the neutral medium {p exrav model in XSPEC, see 
Magdziarz fc Zdziarski 1995 ) and an intrinsically narrow 
emission line at the energy 6.4 keV. No ionization effects 
were taken into account. The binary system inclination 
angle was fixed a,t 9 — 45° (e.g. Zdziarski et al. 199q ), 
the iron abundance - at the solar value. In such a model 
the amplitude of the reflected component is characterized 
by the reflection scaling factor R, which is an approxi- 
mate measure of the solid angle subtended by the reflector, 
R ~ 51/27r. In the simplest geometry of an isotropic point 
source above the infinite reflecting plane, the reflection 
scaling factor R is equal to 1. In order to approximately 
account for smearing of the reflection features due to e.g 
relativistic and ionization effects the reflected continuum 
and the fluorescent line were convolved with a Gaussian. 
Its width was a free parameter of the fit. The uncertain- 
ties in Table |l| represents la confidence intervals for the 
model parameters. The error bars on the values of equiv- 



alent width of the line were calculated by the propagation 
of errors from line flux value. 

The power spectra of GX339-4 in the low spectral state 
feature a prominent QPO peak which frequency varies 
typically between ^^ 0.1 and ~ 0.5 Hz (Fig.n^). We there- 
fore used its centroid to parameterize the characteristic 
noise frequency. The power spectra were approximated 
with a model consisting of two band limited noise com- 
ponents (Lorentzians, centered at zero frequency) and the 
comparably narrow Lorentzian profile (QPO peak). 

The Fourier frequency resolved spectra were obtained 
following the prescription of Revnivtsev, Gilfanov & 
Churazov, 1999 and were approximated with the same 
model as averaged spectra except that the width of the 
Gaussian used to model the smearing of the reflection fea- 
tures was fixed at the value of 0.7 keV. 

3. Results and their uncertainties 

The results of the energy and power spectra approxima- 
tion are presented in Table [^ and Fig.|[ As can be seen 
from Fig.0 the main temporal and spectral parameters 
- characteristic noise frequency, slope of the underlying 
power power law and amplitude of the reflected compo- 
nent ~ change in a correlated way. A steepening of the 
spectrum is accompanied with an increase of the refiection 
amplitude and an increase of the QPO centroid frequency. 
Such a behavior is very similar to that found in Paper I 
for GygX-1. 

The spectral model is obviously oversimplified and 
does not include several important effects such as ion- 
ization of the reflecting media, deviations of the primary 
emission spectrum from the power law, exact shape of the 
relativistic smearing of the reflection features etc. These 
effects might affect the best fit parameters and could 
lead to appearance of artificial correlations between them. 
Particulary sensitive to the choice of the spectral model is 
the refiection scaling factor R ~ r2/27r. As is well known 
there is some degeneracy between the amplitude of refiec- 
tion R and the photon index F of the underlying power 
law determined from the spectral fits, especially if the 
spectral fitting was done in a limited energy range (e.g. 
^dziarski et al. 199£ ). This degeneracy might result in a 
slight positive correlation between the best fit values of 
R and F which is in part due to statistical noise and, in 
part, due to inadequate choice of the spectral model. The 
statistical part of this degeneracy is illustrated in Fig.S by 
a 2-dimensional confidence contour for one of the points 
in the R — T plane. As can be easily seen from Fig.0 it 
is correctly represented by the error bars assigned to the 
points. 

In order to estimate contribution of the second, sys- 
tematic part of the i? — F degeneracy we compare two 
pairs of observations with different and close best fit val- 
ues of reflection factors R. We plot in Fig.|| the ratios of the 
count spectra for each pair. As is clearly seen from Fig|| 
the spectrum having larger best fit value of the refiection 
scaling factor (and the equivalent width of the line) shows 
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Fig. 2. The dependences of reflection scaling factor and the QPO frequency on the slope of the primary power law. 
The filled oval shaped region on the left panel shows a typical 1-cr confidence contour. 

Table 1. The list of observations and parameters of the energy and power density spectra approximation. 



Obs.ID 


Date 


Time, UT 


Exp." 


r 


R ~ n/2-n 


EW, eV 


g", keV 


f QPO, Hz 


X'(40dof) 


10068-05-01-00 


17/10/96 


02:12-05:07 


5550 


1.68 ±0.01 


0.29 ±0.04 


143 ± 30 


0.64 ±0.09 


0.09 ±0.01 


26.5 


10068-05-02-00 


29/10/96 


22:33-00:45 


5530 


1.69 ±0.01 


0.30 ± 0.04 


120 ± 30 


0.62 ±0.14 


0.11 ±0.01 


37.6 


20056-01-01-00 


05/04/97 


08:36-09:15 


2118 


1.79 ±0.01 


0.41 ±0.05 


173 ± 40 


0.85 ±0.11 


0.37 ±0.03 


29.0 


20056-01-02-00 


10/04/97 


11:47-12:28 


2071 


1.81 ±0.01 


0.41 ±0.04 


184 ± 32 


0.75 ±0.08 


0.39 ±0.03 


38.7 


20056-01-03-00 


11/04/97 


13:25-14:06 


2122 


1.81 ±0.01 


0.45 ±0.05 


202 ± 38 


0.83 ±0.09 


0.36 ±0.03 


36.2 


20056-01-04-00 


13/04/97 


20:09-20:50 


2282 


1.82 ±0.01 


0.46 ± 0.05 


213 ±41 


0.92 ±0.09 


0.44 ± 0.02 


32.5 


20056-01-05-00 


15/04/97 


20:41-21:20 


2049 


1.84 ±0.01 


0.47 ±0.05 


218 ±37 


0.89 ±0.08 


0.53 ±0.03 


27.9 


20056-01-06-00 


17/04/97 


23:25-00:01 


2033 


1.84 ±0.01 


0.49 ±0.05 


225 ± 37 


0.92 ±0.07 


0.49 ±0.04 


34.2 


20056-01-07-00 


19/04/97 


22:20-23:08 


2045 


1.86 ±0.01 


0.53 ±0.05 


216 ±33 


0.79 ±0.07 


0.57 ±0.02 


34.2 


20056-01-08-00 


22/04/97 


21:53-22:29 


1988 


1.86 ±0.01 


0.52 ±0.05 


227 ± 38 


0.89 ±0.07 


0.66 ±0.01 


41.7 


20181-01-01-01 


03/02/97 


15:56-19:09 


6622 


1.80 ±0.01 


0.46 ±0.04 


191 ± 30 


0.84 ±0.07 


0.34 ±0.02 


37.7 


20181-01-01-00 


03/02/97 


22:26-01:17 


5452 


1.79 ±0.01 


0.45 ± 0.04 


196 ± 30 


0.83 ±0.07 


0.35 ±0.02 


40.6 


20181-01-02-00 


10/02/97 


15:49-20:22 


10535 


1.79 ±0.01 


0.48 ±0.03 


190 ± 29 


0.86 ±0.07 


0.32 ±0.01 


22.8 


20181-01-03-00 


17/02/97 


18:28-23:46 


11412 


1.79 ±0.01 


0.45 ± 0.03 


176 ± 26 


0.79 ±0.07 


0.32 ±0.01 


24.3 


20183-01-01-01 


08/02/97 


14:20-20:25 


13702 


1.80 ±0.01 


0.47 ±0.03 


185 ± 26 


0.80 ±0.06 


0.37 ±0.02 


34.3 


20183-01-02-00 


14/02/97 


00:18-06:33 


9780 


1.78 ±0.01 


0.40 ±0.03 


177 ± 24 


0.75 ±0.06 


0.34 ±0.01 


21.4 


20183-01-02-01 


14/02/97 


14:20-21:22 


5779 


1.78 ±0.01 


0.43 ±0.04 


182 ± 30 


0.81 ±0.08 


0.30 ±0.02 


30.0 


20183-01-03-00 


22/10/97 


03:00-05:52 


6556 


1.84 ±0.01 


0.51 ±0.04 


203 ± 35 


0.90 ±0.08 


0.42 ± 0.02 


38.5 


20183-01-04-00 


25/10/97 


03:22-06:00 


5385 


1.79 ±0.01 


0.42 ± 0.04 


160 ± 28 


0.72 ±0.09 


0.39 ±0.01 


35.8 


20183-01-05-00 


28/10/97 


18:08-22:13 


6534 


1.76 ±0.01 


0.37 ±0.04 


151 ±27 


0.72 ±0.09 


0.36 ±0.01 


28.2 


20183-01-06-00 


31/10/97 


19:41-22:10 


4621 


1.74 ±0.01 


0.34 ±0.04 


129 ± 29 


0.63 ±0.12 


0.25 ±0.02 


44.2 


20183-01-07-00 


03/11/97 


20:35-23:48 


7186 


1.74 ±0.01 


0.38 ±0.04 


135 ± 26 


0.67 ±0.10 


0.29 ±0.02 


30.7 



" - Dead time corrected value 

' - The width of the Gaussian used to model the smearing of the reflection features. 



more pronounced reflection signatures - the fluorescence 
line at ~ 6.4 keV followed by the absorption edge and in- 
crease due to the Compton reflected continuum at larger 
energies. Thus we conclude that although the best flt val- 
ues of the model parameters might not represent the exact 
values of the physically important quantities, our spectral 
model does correctly rank the spectra according to the 



amplitude of the reflection signatures and the correlations 
shown in Fig.^ are not artiflcial. 

The Fourier frequency resolved spectra illustrate the 
energy dependence of the amplitude of X-ray flux vari- 
at ion at a given tim e sc ale. As was stressed out by 



Revnivtsev et al. 1999 and Gilfanov et al. 2000t interpre- 



tation of the Fourier frequency resolved spectra in general 
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Fig. 4. The ratio of the frequency resolved spectra to a 
power law model with photon index of a=1.8. The 9-25 
Hz spectrum was multiplied by 0.7 for clarity. 



Fig. 3. Ratios of the observed counts spectra with differ- 
ent best fit values of the reflection scaling factor R and 
photon index F. Solid circles show the ratio of the spec- 
trum averaged over April 13-22, 1997 observations, having 
large best fit value of the reflection {R = 0.5, F = 1.84) 
to the average spectrum of Oct. 17 and 19, 1996 with the 
smallest best flt value of reflection and hardest power law 
photon index {R = 0.3, F = 1.69). The ratio was multi- 
plied by i<;~009 for clarity. For comparison the open circles 
show the ratio of the counts spectra of Oct. 25, 1997 and 
Feb. 14, 1997 having close values of the best fit parameters 
(i? = 0.43, F = 1.78 and R = 0.42, F = 1.79 respectively) 
and separated by the similar period in time as the first 
ratio. 



is not straightforward and requires certain a priory as- 
sumptions to be made. One of the areas where it can 
be efficiently used and which at present can not be ac- 
cessed by the conventional spectroscopy is studying the 
fast variability of the reflected emission. Indeed, varia- 
tion of the parameters in the Comptonization region, for 
instance, lead to variations of the spectral shape of the 
Comptonized radiation which would be imprinted in the 
Fourier frequency resolved spectra. Their shape, however, 
might differ significantly from any of the Comptonized 
spectra they resulted from and no easily interpretable re- 
sults could be obtained via the conventional spectral fits. 
The shape of the reflection signatures, on the other hand, 
and especially that of the fluorescent line, is generally 
subject to significantly less variations. Therefore they can 
be easily identified in the frequency resolved spectra and 
their amplitude can be measured. Absence or presence of 
the reflection signatures in a Fourier frequency resolved 
spectrum would signal absence or presence of variations 
of the reflected emission at the given frequency. Their 
amplitude would in principal measure amplitude of the 
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Fig. 5. The Fourier frequency dependence of the equiva- 
lent width of the fluorescent line and the reflection scaling 
factor R. The upper limit on the reflection scaling factor in 
the 9-25 Hz frequency range equals to 1.5 and not shown 
in the graph. 



variations of the reflected flux relative to the variations 
of the primary emission. With that in mind we show in 
Fig.H the Fourier frequency resolved spectra of GX339-4 
in two frequency ranges. Signiflcant decrease of the am- 
plitude of the reflection features with frequency is appar- 
ent. Quantitative dependence of the reflection amplitude 
on the Fourier frequency is shown in Fig.^. Within the 
available statistical accuracy this dependence is qualita- 
tively and quantitatively similar to that found for Cyg 
X-1 ( iRevnivtsev et al. 1999| ). 
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4. Discussion 

We analyzed 23 observations of GX339-4 with 
RXTE/PCA performed from 1996-1997 during the 
low spectral state of the source. Using simple spectral 
model we found that the pattern of temporal and spectral 
variability of GX339-4 is confirming the previous findings 
of [Ueda et al 1994 obtained with the help of GINGA 
observatory, and it is n early identical to that of Cyg X-1 
( Revnivtsev et al. 1999 , Paper I). This indicates that 
such a pattern might be common for the accreting black 
holes in the low spectral state. In particular: 

1. The characteristic noise frequency, slope of the 
Comptonized spectrum and amplitude of the reflected 
component change in a correlated way. Increase of the 
noise frequency is accompanied by increase of the am- 
plitude of the reflected component and steepening of 
the Comptonized spectrum. 

2. Fourier frequency resolved spectral analysis showed 
that the short term variations of the reflected flux are 
suppressed in comparison with variations of the pri- 
mary Comptonized flux at frequencies above ^1 — 10 
Hz. 

As it was discussed e.g. in Gilfanov et al. 1999, 



Zdziarski et al. 1999 the correlation of the spectral in- 



dex (the photon index T of the underlying power law) 
with the amplitude of the reflected component hints 
on the possible relation between the solid angle, sub- 
tended by the reflector, and the influx of the soft pho- 
tons to the Comptonization region. Existence of such 
relation suggests that the reflecting medium is the pri- 
mary source of the soft seed photons to the Comptonizing 
region. It could be explained e.g. within the frame- 
work of t he disk-spheroid mod els of the accretion flow 
(see e.g. Poutanen et al. 1997 ) (a hot quasi-spherical 
Comptonization region near the compact object sur- 
rounded by an optically thick cold accretion disk termi- 
nating at some radius -Kin). In such geometry the decrease 
of the inner radius of the optically thick cold accretion disk 
should result in: 1) the decrease of the temperature of the 
inner V|nt region, because of the inc.rea.se of the influx of 



The frequency resolved spectral analysis was intro- 
duced in Revnivtsev et al. 1999. In that paper we showed 
that the reflected component in the spectrum of Cyg X-1 
is less variable than the underlying continuum at the fre- 
quencies above ~1 Hz. The similar behavior, but with less 
statistical significance was found in the case of GX 339-4. 
Similar to the discussion in Revnivtsev et al. 199E we can 



assume here that the time variations of the reflected com- 
ponent in the spectrum of GX 339-4 could be smeared 
out by the finite light-crossing time of the refiector 



(see more extended discussion in Gilfanov et al. 2000b). 
Alternatively, the observed behavior could be explained 
by non-uniformity within the comptonizing region. For ex- 
ample, if the short time scale variations appear in geomet- 
rically inner part of the accretion flow (hot spheroid) and 
give a rise to significantly weaker, if any, reflected emis- 
sion than the longer time scale events (originating in the 
outer regions) then we will see no reflection features at 
high Fourier frequencies. In turn, the smaller amplitude 
of the reflection in the inner regions of the accretion flow 
might be caused by the screening of the reflector from the 
innermost regions by the outer parts of spheroid. 
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